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Background 


m Overvoltage studies in power systems normally performed 
using EMIP tools (ElectroMagnetic Transients Programs) 
B EMTP-rv 
m PSCAD 
m ATP 


E EMIP tools do not include wide-band, transformer 
models: 


B Accurate for low frequencies 


m Not suitable for simulation of high-frequency transformer- 
network interactions. 
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Scope 
m [erminal equivalent model suitable for high-frequency 


transformer network interaction studies 
B Linear effects only 


m Interfacing the model with EMTP-type tools 


m Input data 
m Frequency sweep measurements, or 
m Manufacturer's white-box model parameters 


m Examples 
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Part I: 
Transformer Terminal Equivalent 
Model 
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Terminal admittance matrix 


«—— terminals 
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Y(o): From measurements or calculations 
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Modeling via rational functions 


I (0) Y,(O) Y (0) = Y, (0) || Vo) 
[,(@) | |n œ) No) Y (o) || V,(@) 
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Vector Fitting 


E Method for fitting a pole-residue model to frequency data 





Y(0)= Y R, +D 


m=1 JO mn a, 
B Robust, efficient, accurate 


m Guaranteed stable poles 
BH Symmetrical model 


m Freely available on internet as Matlab function 


“B. Gustavsen and A. Semlyen,"Rational approximation of frequency domain 
responses by vector fitting", IEEE Trans. Power Delivery, July 1999. 
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Passivity 


m [he model extracted by VF is not necessarily passive 


m Loss of passivity often leads to unstable time domain 
simulations 
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Passivity enforcement y, y ke -+p 


m=l JO — da, 


m Perturb model's parameters so that the model becomes 
passive" 
B Perturb elements in {R,,},D 
B make change to Y(s) minimal 


N 
AY =Y An +AD=0 








m=l JØ anm 
N AR 
eig(Re{Y+ > ——1)>0 L | 
mai JO—a, > A, AX =0 Constrained LS 
eig(D+ AD) > 0 B, Ax «c (e.g. QP) 


Ax: elements of {AR}, AD 


"B. Gustavsen, and A. Semlyen, “Enforcing passivity for admittance matrices approximated 
by rational functions’, IEEE Trans. Power Systems, vol. 16, no. 1, pp. 97-104, Feb. 2007. 
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Part Il: 
Interfacing the model with EMTP tools 
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Including model in EMTP type program 


N 
Y(o)=> Rn DHE 


m=l JOU, 





Lumped circuit 


equivalent 
EMTP time domain 


or simulation 


Norton equivalent 
w/ controlled source 
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1. Lumped circuit equivalent 





med 
| 

uM ja AN 
mel 


laila — Xp 


EMTP branch cards 
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2. Norton equivalent w/ controlled source 








yo == y(t)=) re” 


Jo-—a 





m Recursive convolution in time domain (fast!) 
X =UX ,+v, 


n 


i(t) = y(t)*v(f) mu» 
Trapezoida 
integration 


… T 
| LC X TYV, 
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EMTP compatibility 





RLC network man 
Norton equivalent 


* User-defined component 


** Convert to state-space block 
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Part Ill: 
Model input data 
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Samples for the model identification 


+ D-+sE 


Y(o) = Yo 





m Frequency sweep measurements, or 


m Calculations from manufacturer's white-box model 
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Calculating Y(o) using frequency sweep 
measurements 
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Columnwise measurement of 
admittance matrix 


Example: 3-terminal device 


Vi l, 
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Frequency sweep measurements 


Network analyzer 
Connection box 
Coaxial cables 
Current sensor 
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Measurement setup (cont.) 


Measuring Y, =i/v 


C — MM WM " pow mms MW : 
T Er. Source Ref. Input | 





Current 
sensor 


Jumper to Å ` 
ground 
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hl Yu Yo 
L, Y, Y, 
la — Yo, Yo 
Ll [Fa fo 
l; Ys, Ys 
A Ya Yo 
Elements of Y(o) 
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ON 
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Measurement of Y (5 Hz — 10 MHz) 


300 kVA transformer 





Vi 
V» 
V3 1 
V, 2 
Vs 3 
Ve 
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Voltage ratio 





Vector Fitting + Passivity Enf. 


100 


W +D-+sE 
Y moder 6 ) = DE a, 





m Voltage ratios 
(computed from Y gata and Y model) 


high-low low- high 


Voltage ratio 


Frequency [Hz] 
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Time domain measurement & simulation 


m Simulation method: Norton equivalent 






A. 
l | 
| | 

peces 
Time [us] 
-1 | 

0 1 2 3 4 5 pa i 
Time [ps] 
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Laj RJ = 








SINTEF 
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100 


150 
Time [us] 


200 
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Alternative approach: 


Calculating Y(o) from manufacturer's 
white-box model paramretrs 
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Winding discretization 


Example =. 














Fictitious transformer example used in a 
CIGRE WG A2/C4.39 (courtesy: Alvaro Portillo) ME 

" x XZ 

F? 35 X— 

LV : se 

1 i = 

100 MVA unit, 230/69 kV, BIL=1050 kV B 
mens — iq os Kid AJN 10.223 HV Winding — “Inzerleaved Disk Type” s = 
c uu T "Continuous Disk Type" win DES The high voltage winding Is a conventional "Interieaved Disk Type" (Figure 10.4). å — 
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White Box Model description 


Lumped circuit parameter model, 
based on spatial discretzation of windings 


Matrices: R, L, G, C 


Basic equations ^ 
e= Tv L: | 
Mutual coupling beween 
# Ta u " 
i, = Ti all inductive branches 


I. +1, +l =I, 
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Calculate terminal admittance matrix 
W.r.t. external nodes 


1. Calculate admittance matrix Y(s) to all nodes, at discrete 
freq. 
Y(o,) 2[(G jo, C) - T' (R+ jo L) T].k-1..K 





2. Introduce groundings and short circuits 


3. Heduce model with respect to external terminals 


You (JO, ) 
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Admittance matrix w.r.t. terminals 47, 70 


Admittance [S] 





100 200 300 400 500 600 700 800 900 
Frequency [kHz] 
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1000 
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Modeling with rational functions 





un (0) = y ym 


+D Rational funcion approximation 


ma JO-a 


N-80 


Admittance [S] 





100 200 300 400 500 600 700 800 900 1000 
Frequency [kHz] 


The information about white-box parameters is hidden ! 
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Internal voltages 


V oxi W Li, Z, Lu 
Via i Le Zp - 
Calculate voltage transfer to internal node as 


LAP (a) — H inter (W) Ved (W) 


H ster (W) — Z C (W) Y... (@ ) 
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Voltage transfer from external nodes 47, 70 to e 
Internal nodes 24 47 58 59 68 69 72 73 82 83 a 


Voltage transfer from node #47 to winding nodes 4 

8 40 

39 

7 38 

37 

36 

6 35 

34 

a 9 33 

a 32 

o 4 31 

CU 30 
= 

3 3 29 





2 26 
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0 200 400 600 800 1000 e 


Voltage transfer from node #70 to winding nodes " 


Voltage [p.u.] 
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Modeling with rational functions 


h;(s) = yt rm „i=1,2 Rational funcion approximation 








m=1 ST (Vector Fitting) 


Voltage transfer from node #47 and #70 to all winding nodes 


N=80 
(per column) 





Voltage [p.u.] 





100 200 300 400 500 600 700 800 900 1000 
Frequency [kHz] 
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Time domain simulation : 






u Voltage at internal nodes : 
t 24 47 58 59 68 69 72 73 82 83 e 
LV —— 971 1.2/50 US : 


Voltage [p.u.] 





o 20 40 60 80 100 
Time [us] 


- N WW A MN O Nê os 
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Simulation in PSCAD 


SINTEF 
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Network interaction study 
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Conclusions 


m [ransformer high-frequency behavior can be 
characterized by its terminal admittance matrix, Y (c) 


m Y(o) can be fitted with a pole-residue model to obtain an 
EMIP-compatible model 
B Lumped circuit equivalent, or 
m Norton equivalent (recursive convolution) 


E Y(o) can be obtained from 
B Measurements, or 
m Manufacturer's white box model parameters 
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